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ABSTRACT 

In the limit of heavy quark mass, the production cross section and polarization of 
quarkonia can be calculated in perturbative QCD. We study the p_i_-averaged pro- 
duction of charmonium states in ttN collisions at fixed target energies. The data on 
the relative production rates of J /if) and xj is found to disagree with leading twist 
QCD. The polarization of the J /if) indicates that the discrepancy is not due to poorly 
known parton distributions nor to the size of higher order effects (if-factors). Rather, 
the disagreement suggests important higher twist corrections, as has been surmised 
earlier from the nuclear target A-dependence of the production cross section. 
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1 Introduction 



Quarkonium bound states formed by heavy quark-antiquark pairs are small nonrela- 
tivistic systems, whose production and decay properties are expected to be governed 
by perturbative QCD. The extensive data available on the inclusive decays of many 
charmonium (cc) and bottomonium (bb) states has been compared with detailed per- 
turbative calculations. The overall agreement between theory and experiment is rea- 
sonably good, taking into account the moderate mass scale |], @, |3|, |J. This work 
has led to self-consistent values of the size of the quarkonium wave function near the 
origin. 

Assuming, then, that we understand the decay of the quarkonium states to per- 
turbative gluon and light-quark final states, we can turn the reaction around and 
consider the photoproduction and hadroproduction of quarkonia. Thus quarkonium 
production becomes a probe of the production mechanism of color-singlet heavy quark 
pairs. This is analogous to the dynamics of lepton pair hadroproduction, where the 
main production mechanism (at lowest order) has been identified as the Drell-Yan 
hard fusion subprocess qq —>■ 7*. Quarkonium production can offer new insights 
into gluon fusion mechanisms; for example, the J/ip and Xi couple to states with 
more than two light partons, such as ggg or qqg. At leading twist, i.e., to leading 
order in I / rug, quarkonium production proceeds through the collision of only two 
partons, one from the projectile and one from the target. Hence an extra gluon or 
quark must be emitted in the leading twist production of J/ip and xi- However, at 
large values of the quarkonium momentum fraction xp, it becomes advantageous for 
two or more collinear partons from either the projectile or target to participate in 
the reaction. Such processes are higher twist, since their rate is suppressed relative 
to ordinary fusion reactions by powers of Aqcd/t^q where Aqcd is the character- 
istic transverse momentum in the incident hadron wavefunction. Nevertheless, de- 
spite the extra powers of I/tuq, the multiparton processes can become dominant at 
(1 — xp) < 0(AQ CD /mQ) since they are efficient in converting the incident hadron 
momentum into high xp quarkonia ||. 
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In leading twist QCD the production of the J/ip at low transverse momentum 
occurs both 'directly' from the gluon fusion subprocess gg — > J/ip + g and indirectly 
via the production of xi an d X2 states^. These states have sizable decay branching 
fractions Xi,2 J/ip + "f of 27% and 13%, respectively. In spite of its relatively small 
branching ratio, the X2 state is expected to give an important contribution to the 
total yield of J/ip's at leading twist, since gg — ► X2 is of lower order in a s compared 
to the competing processes. Early comparisons |7|, §, |Hj with the total J/ip cross 
section data indicated rough agreement with the model predictions. Nevertheless, 
the cross sections for direct J/ip and xi production were predicted [jy]] to be too low 
compared to the data (TJ, |TJ, [TJ, |T5[ 0. 



More recent E705 and E672 data |17|, |18| on the production fractions of the various 
charmonium states have confirmed that there is a clear discrepancy with the leading 
twist QCD prediction. The leading twist calculations which we present in this paper 
show that the predicted ratio of direct J/ip production in ttN collisions compared 
to the X2 production is too low by a factor of about 3. In addition, the ratio of 
Xi production to X2 production is too low by a factor of 10. A similar conclusion 
has been reached in ||, where possible explanations in terms of uncertainties in 
the partonic cross sections (very different fT-factors for the various processes) or 
unconventional pion parton distributions are discussed. Less data is available for 
proton-induced charmonium production, but a discrepancy between leading twist 
QCD and experiment appears likely also in that case. 



The wealth of data from the NA3 experiment at CERN |]19[ and the Chicago-Iowa- 
Princeton [|2D| and E537 experiments |^TJ at FermiLab on the angular distribution of 



the muons in the decay J/ip — > fi + fi provides an even more sensitive discriminant of 



different production mechanisms |22L G3L E4L RR 26, 27, E 



The polarization of the cc, and hence that of the charmonium bound state [23|. 

can at leading twist be calculated from perturbative QCD. Furthermore, in the heavy 

quark limit, the radiative transition xj —> J/V + 7 preserves the quark spins, i.e., it is 

an electric dipole transition. Hence the polarization also of indirectly produced J/^'s 

^^At high transverse momentum, one also has to take into account production through quark and 
gluon fragmentation ||. 
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can be calculated. We find that even if the relative production rates of the J/if), Xi 
and X2 are adjusted (using X-factors) to agree with the data, the J/if> polarization 
data is still not reproduced. 

We shall argue that a possible explanation for the underestimate of the J /if) and 
Xi cross sections is that more than one parton from either the projectile or target 
participates in the collision, so that no additional gluon needs to be emitted. Similar 
higher twist effects are known to become important at high xp in lepton pair pro- 



duction |29|, |3(], [31], [32]]. In Ref. [|3(J it is shown that higher twist contributions can 
explain the large azimuthal cos0 and cos 20 correlations seen in the ttN — > yU + /i~ 
data. There are also previous indications from the non-factorizing anomalous nuclear 



that higher twist effects 



target dependence of the J / if) cross section |33], [19], |34], [35], 
are considerably larger in J /if) production than in lepton pair production, and that 
they persist down to low xp. 



2 Production rates of ip and \J states at leading 
twist 

In this section we calculate J /if) production in 7rN interactions at leading twist and 
to lowest ordeiQ in a s . Higher order corrections in a s and relativistic corrections to 
the charmonium bound states are unlikely to change our qualitative conclusions at 
moderate Xp. Contributions from direct J /if) production, as well as from indirect 
production via Xi an d Xi decays, are included. Due to the small branching fraction 
Xo — ¥ J/4> + 7 °f 0.7%, the contribution from xo to J /if) production is expected (and 
observed) to be negligible. Decays from the radially excited 2 3 Si state, if)' — > J/if)+X, 
contribute to the total J /if) rate at the few per cent level and also will be ignored 
here. 

Since the if)' is formed directly, its production allows an important cross check on 
the use of charmonium states to study the production mechanism. At high energies, 
the charmonium bound state forms long after the production of the compact cc pair 
(the formation time Tf orm ~ 2E\ & T b / AM 2 ). Thus the ratio of if)' to direct J /if) produc- 
2 Thus we do not include subprocesses like qg — * \2l (which is subleading to gg — > xz)- 
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tion can depend only on the relative magnitude of their wave functions at the origin. 
More precisely (see, e.g., 0]), 

a(ip') T(f -> e+e" 



e T e" 



0.24 ±0.03 



CT dk (J/lp) T(J/^ 

where cr^ r {J/i})) is the cross section for direct production of the J/ip. The ratio 
([I]) should hold for all beams and targets, independent of the size of the higher twist 
corrections in producing the pointlike cc state. The energy should be large enough for 
the bound state to form outside the target. The available data is indeed compatible 
with (|J). In particular, the E705 value |17| is about 0.24 (see Table 1). The anomalous 
nuclear target A-dependence observed for the J/ip is also seen for the ifi' so that 
the ratio (]!]) is indeed independent of A. 





a(ip') [nb] 


CTdir^/V^rib] 




7T+ 


22 ±5 


97 ± 14 


0.23 ±0.07 




25 ±4 


102 ± 14 


0.25 ±0.05 


V 


20 ±3 


89 ± 12 


0.23 ±0.05 



Table 1: Production cross sections for ■*/>', direct J/ip and their ratio in ir + N, tt N 



and pN collisions. The data are from Ref. [17 



The ttN — > X2 + X production cross section to lowest order and twist is 

a{7fN -> X2 + X; x F > 0) = r / F g/n (xi)F g/N (T/x l )a (gg -> xz) (2) 

J t/t X\ 

where r = M^Js and the quantity a (gg -> xs) = l ^ 2 a 2 s \B! P (ti)\ 2 /M^ E3- We re " 
strict the %2 momentum range to the forward CM hemisphere (x F > 0) in accordance 
with the available data. 

The direct ttN — > J/ip + X cross section is similarly given by 

a(7riV — >• J/^ ± X; x F > 0) = / da;i / dx 2 L dtF g /^(x 1 )F g / N (x 2 ) 

Jt Jt/x 1 Jt min 

x^(gg^ J/^ + g) (3) 
at 



where t is the invariant momentum transfer in the subprocess, and 



'x 2 M 2 j 



max 



Xi + x 2 
5 



(4) 



Eq. (^) also applies to the tcN — > \i + X reaction, in which case a sum over the 
relevant subprocesses gg — > Xi9, 9Q ~ > XiQi 9Q ~> XiQ an d — ► Xi<? is necessary. The 



differential cross sections da/dt for all subprocesses are given in [10, 57]. In Table 
2 we compare the Xi production cross section and the relative rates of direct J/ip 
and xi production at Ei a b = 300 GeV with the data of E705 and WAll on n^N 
collisions at Eiab = 300 GeV and 185 GeV 0, [L7|]. We use the parton distributions 
of Ref. [^, ^] evaluated at Q 2 = M 2 , where M is the mass of the charmonium 
state in question. We take a s = 0.26 for all states and use |.Rs(0)| 2 = 0.7 GeV 3 , 
\R' P (0)/M\ 2 = 0.006 GeV 3 0. 





<*{X2) [nb] 


^dir(</AA)MX2) 


<Kxi)Mxa) 


Experiment 


188 ±30 ±21 


0.54 ±0.11 ±0.10 


0.70 ±0.15 ±0.12 


Theory 


78 


0.17 


0.067 



Table 2: Production cross sections for xi, % 2 and directly produced J/^ in n N 



collisions at 300 GeV. The data from Ref. [|14| , p!7|| include measurements at 185 and 
300 GeV. 

The X2 production rate in QCD agrees with the data within a 'ZT-factor' of or- 
der 2 to 3. This is within the theoretical uncertainties arising from the J ftp and \ 
wavef unctions, higher order corrections, parton distributions, and the renormaliza- 
tion scale. A similar factor is found between the lowest-order QCD calculation and 
the data on lepton pair production PT| , [4"2"fl . On the other hand, Table 2 shows a 
considerable discrepancy between the calculated and measured relative production 
rates of direct J/ip and xi compared to xi production. A priori we would expect 
the i^-factors to be roughly similar for all three processes. It should be noted that 
there is a kinematic region in the J/ip and Xi processes where the emitted parton 
is soft (in the rest frame of the charmonium), and where perturbation theory could 
fail. However, the contribution from this region is numerically not important (there 
is actually no infrared divergence). Hence one cannot hope to boost the cross section 
significantly by multiplying the soft parton contribution by any reasonable factor. 
Moreover, the same soft parton region exists in charmonium decays, where analogous 
disagreements with data are absent. It should also be noted that the contribution to 
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Xi production from the qq — > Xi + 9 subprocess is apparently singular at s — M 2 ^ 



due to a breakdown of the non-relativistic approximation of the bound state fTT |. 
The divergence is cancelled once one takes into account higher Fock states In 
agreement with Ref. |33| we find that the cross section is insensitive to the value 



of the cutoff parameter excluding the soft gluon region. 

We conclude that leading twist QCD appears to be in conflict with the observed 
rate of direct J/ip and Xi production. Although in Table 2 we only compared our 
calculation with the E705 and WAll n~N data, this comparison is representative of 
the overall situation (for a recent comprehensive review see ||). 



3 Polarization of the J/ip 



The polarization of the J/ip is determined by the angular distribution of its decay 
muons in the J/ip rest frame. By rotational symmetry and parity, the angular distri- 
bution of massless muons, integrated over the azimuthal angle, has the form 

da 



dcosO 



oc 1 + A cos 2 6 (5) 



where we take 9 to be the angle between the fi + and the projectile direction (i.e., 
we use the Gottfried-Jackson frame). The parameter A can be calculated from the 
cc production amplitude and the electric dipole approximation of radiative x decays. 
Earlier calculations of the polarization in hadropro duct ion p2 , [24| , |25|j were based on 



general effective couplings of the quarkonia and partons rather than the perturbative- 
QCD matrix elements which we shall use. 

The electric dipole approximation of the radiative decay xj ~* ipl is exact in the 
heavy quark limit, i.e., when terms of 0(E 1 /m ( ^) are neglected. As a consequence, the 
heavy quark spins are conserved in the decay, while the orbital angular momentum 
changes. This spin conservation may also be derived from Heavy Quark Symmetry 



44]. The validity of the electric dipole approximation for xj radiative decays has 



been verified experimentally [[4l| . 
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The amplitude for the \2 production subprocess g(/ii)g(/i 2 ) — ► cc —>■ X2{J Z ) is, 
following the notation of Ref. |4(| . 



A{J Z = ±2) = da^Oj^e^ 



x [1 =p (^2 — pt>i) cos $ — yUi/i 2 cos 2 i? - 5 MlM2 sin 2 1?] (6) 
A(J Z = ±1) = Aa a H P {Q) s l^wy& 

x [//i - /i 2 T 2/ii// 2 cos i? ± 25 MlAt2 cos #] (7) 

A(J z = 0) = 4a s R P {0)^pV6 6^ 2 sin 2 $ (8) 

where are the polar and azimuthal angles of the beam gluon in the Gottfried- 
Jackson frame; d = if the transverse momenta of the incoming gluons are neglected. 
In this expected for physical, transversely polarized gluons with /ii j2 = ±1, 

the amplitude for \2 production with J z — jj,x — \ii — ±1 vanishes. Surprisingly, the 
amplitude for J z = also vanishes when $ = 0. Hence the \2 is at lowest order 
produced only with J z = ±2. In this polarization state the spin and orbital angular 
momenta of its constituent charm quarks are aligned, S z = L z = ±1. Since S z is 
conserved in the radiative decay % 2 — > J/tj) + 7, it follows that J Z (J /ip) = S z = ±1 
(L = for the J ftp). Thus the J/ip's produced via X2 decay are transversely polarized, 
i.e., A = 1 in the angular distribution ([5]). This result is exact if both the photon 
recoil and the intrinsic transverse momenta of the incoming partons are neglected. 
Smearing of the beam parton's transverse momentum distribution by a Gaussian 
function exp [— (fcj_/500 MeV) 2 ] would reduce A to ~ 0.85. 

From the gg — > J/ip + g amplitude we find for direct J/ip production, ttN — > 
J/ip + X -> + x > 

1 da 3 r dx\d%2 



x 



D i;li r/.r, (i \n J (xx + X2) i F 9/^ F 9/^ 

0u + 0oo + (Qii ~ Qoo) cos 2 9 (9) 



where is the J/ip — > branching fraction, xp = 2p^/y/s is the longitudinal- 

momentum fraction of the J/ip, and 9 is the muon decay angle of Eq. @. The 
density matrix elements Qn, g 00 are given in the Appendix. 
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For the ttN + 7 + + 7 + X production process we 

obtain similarly 

1 da 3 . 1 \\r^ f dx\dxi 



-Br(xi -> V7) E / T^rflf ^(^O^j/JV^) 



B^dxpd cos 9 1287T (xi + 



x 



0oo + 3^n + (&>o - 0ii ) cos , (10) 



where the density matrix elements for ij = gg, gq, gq and qq scattering are again 
given in the Appendix. 

In Fig. la we show the predicted values of the parameter A of Eq. (|5]) in the 
Gottfried-Jackson frame as a function of xp, for the direct J/ip and the Xi,2 J/ip+j 
processes separately. Direct J /if) production gives A ~ 0.25 in the moderate xf region, 
whereas production via xi results in A ~ —0.15. The dashed lines indicate the effect 
of a Gaussian smearing in the transverse momentum of the beam partons. 

The A(xi?)-distribution obtained when both the direct and indirect J /if} production 
processes are taken into account is shown in Fig. lb and compared with the Chicago- 
Iowa-Princeton [2(J and E537 |H data. Our QCD calculation gives A ~ 0.5 for 



xp < 0.6, significantly different from the measured value A ~ 0. The E537 data gives 
A = 0.028 ± 0.004 for xf > 0, to be compared with our calculated value A = 0.50 in 
the same range. 

The discrepancies between the calculated and measured values of A are one further 
indication that the standard leading twist processes considered here are not adequate 
for explaining charmonium production. The J/if> polarization is particularly sensi- 
tive to the production mechanisms and allows us to make further conclusions on the 
origin of the disagreements, including the above discrepancies in the relative produc- 
tion cross sections of J/ip, Xi an d X%- If these discrepancies arise from an incorrect 
relative normalization of the various subprocess contributions (e.g., due to higher 
order effects), then we would expect the J/ip polarization to agree with data when 
the relative rates of the subprocesses are adjusted according to the measured cross 
sections of direct J /if), Xi an d Xi production^. The lower curve in Fig. lb shows the 

3 In the case of Drell-Yan virtual photon production, it is known that higher order corrections 
do not change the 7* polarization significantly | |f7| , which makes it plausible to represent these 
corrections by a simple multiplicative factor that does not affect the polarization of the photon. 
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effect of multiplying the partial J / ip cross sections with the required X-factors. The 
smearing effect is insignificant as shown by the dashed curve. The A parameter is still 
predicted incorrectly over most of the xf range. 

A similar conclusion is reached (within somewhat larger experimental errors) if 
we compare our calculated value for the polarization of direct J / ip production, shown 
in Fig. la, with the measured value of A for ip' production. In analogy to Eq. 
dU), the ip' polarization data should agree with the polarization of directly produced 
J/ift's, regardless of the production mechanism. Based on the angular distribution 
of the muons from ip' — > decays in 253 GeV ir~W collisions, Ref. |32] quotes 

= 0.02 ± 0.14 for xp > 0.25, appreciably lower than our QCD values for direct 
J/^'s shown in Fig. la. 



4 Discussion 

We have seen that the J/ip and xi hadroproduction cross sections in leading twist 
QCD are at considerable variance with the data, while the \2 cross section agrees with 
measurements within a reasonable i^-factor of 2 to 3. On the other hand, the inclusive 
decays of the charmonium states based on the minimal perturbative final states (gg, 
ggg and qqg) have been studied in detail using perturbation theory [p], ||, El H> anc ^ 
appear to be fairly well understood. It is therefore improbable that the treatment of 
the cc binding should require large corrections. This conclusion is supported by the 
fact that the relative rate of ij}' and direct J/ip production (Eq. ([]])), which at high 
energies should be independent of the production mechanism, is in agreement with 
experiment. 

In a leading twist description, an incorrect normalization of the charmonium pro- 
duction cross sections can arise from large higher order corrections or uncertainties in 
the parton distributions || . Even if the normalization is wrong by as much as a factor 
of 10, such a i^-factor would not explain the J ftp polarization data. Thus a more 
likely explanation of the discrepancy may be that there are important higher-twist 
contributions to the production of the J/ip, if >> and Xi- 

The direct J/ip and xi subprocesses require, at leading order and twist, the emis- 
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sion of a quark or gluon, e.g., gg — ► J/ip + g. This implies a higher subenergy \/d 
for these processes compared to that for the \2i which can be produced through sim- 
ple gluon fusion, gg — > \2- It is then plausible that a higher twist component which 
avoids the necessity for gluon emission is more significant for the J/ip and the Xi than 
it is for the \2 (or for lepton pair production, qq — > 7*). If either the projectile or the 
target contributes two partons rather than one, no emission of a parton is required 
H] : {99) + 9 ~^ J/ip- Similar production mechanisms are considered in Ref. 0. 

Taking two partons from the same hadron is a higher twist process, and as such is 
suppressed by a factor of 0(k 2 Q CD /m 2 c ). This factor describes the probability that the 
two partons are within a transverse distance of 0(l/m c ), as required if both of them 
are to couple to the same cc pair. For the J /if) and Xi, this suppression is compensated 
by the fact that the subprocess energy s can be equal to the charmonium mass since 
no parton needs to be emitted. Finding two softer gluons in a hadron may also be 
more probable than the probability for one gluon carrying the full momentum. 

In the xp — > 1 limit, important higher twist effects are expected |29|, [3(| and 
observed Ej| also in the muon pair production process, ttN — ► fi + ^ + X. In effect, 
both valence quarks in the pion projectile must be involved in the reaction if the full 
momentum is to be delivered to the muons. The higher twist effect manifests itself in 
the angular distribution of the muons: the polarization of the virtual photon changes 
from transverse to longitudinal at large xp. Thus the photon tends to carry the same 
helicity as the pion in the xp — > 1 limit. It is natural to expect the higher twist effects 
to be similarly enhanced in J/ip production at large Xp. As seen in Fig. lb, the data 
does indeed show a remarkable turnover in the polarization of the J /if) for xp ^ 0.8, 
with the fastest J/^'s being longitudinally polarized. In contrast to the lepton pair 
production case, the evidence for higher twist effects persists, as we have seen, for 
J/if>'s produced even at lower momentum fractions. 

It has recently been pointed out ]|9], that there is also a large discrepancy 
between the leading twist QCD prediction and data for large p± charmonium pro- 
duction. At leading twist in l/p± and in 1/m 2 . the dominant source of "prompt" 
J/^'s (i.e., those not due to B — > J/if> + X decays) is predicted to be the radiative 
decays of P-wave charmonia, \ c — > J /if} + 7. The J/if> cross section obtained this 
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way is consistent with the data within a factor ~ 2. The actual Xc production cross 
section has not yet been measured. However, the prediction for tp' production (which 
cannot be produced via radiative decays) is too low by a factor ~ 30 when compared 
to the data. At the same time, the experimental ratio of the tp' and total (prompt) 
J ftp cross section is consistent with the universal ratio of Eq. (|l|). This suggests 
that a major part of the prompt high p± J ftp's are produced directly, rather than via 
X c — ► J/ip+ry decays. Since the shape of the ^-distribution of the measured J ftp and 
tp' cross sections is in agreement with the leading twist prediction ||9|, |50]], the large 
higher twist corrections are likely to reside in the g —>■ J ftp and c — > J ftp fragmenta- 
tion vertices, and thus be of 0(l/m^) rather than of 0(l/p 2 ± ). This is consistent with 
our conclusions based on the low p± charmonium data. The much larger discrepancy 
at high p± is qualitatively expected since the high z region of the fragmentation is 
emphasized due to the "trigger bias" effect. As discussed above, the higher twist 
"intrinsic charm" mechanism is particularly important at high momentum fraction in 
either the projectile or fragmenting parton systems. 

Additional independent evidence for higher twist effects in J / tp production is also 
reflected in the nuclear target A-dependence of the cross section. In lepton pair 
production, the cross section is very closely linearly dependent on A (apart from a 
small deviation at the largest xp [34|). J ftp production, on the other hand, shows 
a nuclear suppression over the whole xp range p6| |. The suppression depends on xp 



rather than on x%, and it is thus possible to conclude that QCD factorization 
must be broken, implying that the effect is due to higher twist terms. 

Further theoretical work is needed to establish that the data on direct J ftp and 
Xi production indeed can be described using a higher twist mechanism of the type 
discussed here. Experimentally, it is important to check whether the J ftp's produced 
indirectly via xi decay are transversely polarized. This would show that X2 produc- 
tion is dominantly leading twist, as we have argued. Better data on real or virtual 
photoproduction of the individual charmonium states would also add important in- 
formation. So far, little is known about the relative size of direct and indirect J ftp 



photoproduction, and the polarization measurements [^T|, [52| are too inaccurate to 
test theoretical predictions |23L |26|, |27|, |28| . 
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In photoproduction one expects less higher twist effects associated with the pro- 
jectile, since in the case of direct photon interactions only a single beam parton (the 
photon itself) is available. However, the target hadron can contribute two gluons. In 
the special case of diffractive J/ip photoproduction pBH , this is in fact expected to be 
the dominant reaction mechanism W< 
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A Density matrix elements 

The density matrix elements for J / ip production in the gg — > J/ip + g subprocess are 
defined as 

Qw! = 256 ^ A (9i92 -> J/i/>(n) + 93)A*{g 1 g 2 -> J/i>(p') + g 3 ), (H) 

^1/^2/^3 

ClC 2 C 3 

where are the helicity and color of the gluon i, and the factor of 1/256 comes 
from averaging over the initial helicities and colors. The diagonal matrix elements 
are found to be 

407r 2 a 3 s \R s (0)\ 2 M 
011 ~ 9[(s - M 2 )(t - M 2 )(u - M 2 )] 2 

x {s 2 (s - M 2 ) 2 + t 2 (t - M 2 ) 2 + u 2 {u - M 2 ) 2 

- 2M 2 [k 2 1± (s 2 + t 2 ) + 2k 1± ■ k 2± s 2 + k 2 2± (s 2 + u 2 )]} , (12) 

407r 2 a^g(0)| 2 M 
^ 00 ~ 9[(s - M 2 ){t - M 2 )(u - M 2 )} 2 

x {s 2 (s - M 2 ) 2 + t 2 (t - M 2 ) 2 + u 2 {u - M 2 ) 2 

- AM 2 [k 2 lz {s 2 + t 2 ) + 2k lz k 2z s 2 + k 2 2z {s 2 + u 2 )}) , (13) 
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where M is the J/ip mass, s, t, u are the subprocess invariants (the carets being 
omitted for clarity), and ki^ are the three-momenta of the beam and target partons 
in the Gottfried- Jackson frame. 

In analogous notation, the diagonal density matrix elements for xi production in 
qq, gq and gg scattering are 



000 



Qil = 


(647r) 2 a 3|^(0)| 2 




9M(s - M 2 ) 4 




QQ 

f?00 _ 


(647r) 2 a3|^( )| 2 




9M{s - M 2 ) 4 




QQ 

m = 


-3(647r) 2 a3|^(0) 


I 2 


72M(t - M 2 ) 4 




QQ 

@oo — 


-3(647r) 2 af|^(0) 


I 2 


72M(t - M 2 ) 4 




Qn — 


96^1^ (0)| 2 




M 3 (Q - M 2 P) 4 





tk 2 1± + (s- M 2 )k 1± • k 2± - uk 2 21 _ + tu 
2tk\ z + 2(a - M 2 )k lz k 2z - 2uk\ z + tu 
— sk 2 _,_ + (u — s)k 1± ■ k 2 j_ + su 
—2sk\ z + 2{u — s)k\ z k 2z + su 



xP 2 
1 



M 2 P 2 (M 4 - 4P) - 2Q(M 8 - 5M 4 P - P 2 ) - 15M 2 Q 2 



-n" 
2 tXIO j 



48tt 2 «3|P^(0)P 



Mstw[(s - M 2 )(t - M 2 ){u - M 2 )] 4 

x {s 2 (s - M 2 ) 2 [k lz g{s, t, u) + k 2z g{s, u, t)] 2 

+ u 2 {u- M 2 ) 2 [k lz g{u, t, s) - {k lz + k 2z )g{u, s, t)] 2 

+ t 2 (t - M 2 ) 2 [(k lz + k 2z )g(t, 8 , u) - k 2z g(t, u, s)} 2 



+ AM 2 (k lx k 2y - k 2x k ly ) 2 



x 



s 2 (s - M 2 ) 2 f 2 (s, t, u) + (s^u) + (s^ t)] } . 



(14) 
(15) 
(16) 
(17) 



(19) 



The density matrix elements g q9 for the processes where a beam quark scatters off a 
target gluon are obtained by changing ki *-> k 2 (and consequently t <-> u) in g 9q . The 
matrix elements for the gq and qg scattering processes are the same as for gq and qg, 
respectively. In deriving g" we made use of the subprocess amplitudes given in Ref. 
[3~7[1 . The functions /, g, P and Q of the invariants are defined as 



f(s,t,u) 



(t — u)(st + tu + us — s z 



(20) 
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g(s,t,u) = (s + t)[st(t-s) + su(u-s)+tu(t-u)}, (21) 
P = s t + tu + us, (22) 
Q = stu. (23) 



15 



References 

[1] W. Kwong, J. L. Rosner, and C. Quigg, Ann. Rev. Nucl. Part. Sci. 37, 325 
(1987). 

[2] L. K6pke and N. Wermes, Phys. Rep. 174, 67 (1989). 
[3] G. A. Schuler, CERN-TH.7170/94 (February 1994). 

[4] G. T. Bodwin, E. Braaten, and G. P. Lepage, Phys. Rev. D46, 1914 (1992); 
Fermilab-PUB-94-073-T (June 1994). 

[5] S. J. Brodsky, P. Hoyer, A. H. Mueller, and W.-K. Tang, Nucl. Phys. B369, 519 
(1992). 

[6] E. Braaten and T. C. Yuan, Phys. Rev. Lett. 71, 1673 (1993). 

[7] C. E. Carlson and R. Suaya, Phys. Rev. D18, 760 (1978). 

[8] J. H. Kiihn, Phys. Lett. B89, 385 (1980). 

[9] C.-H. Chang, Nucl. Phys. B172, 425 (1980). 

[10] R. Baier and R. Riickl, Z. Phys. C19, 251 (1983). 

[11] V. Barger and A. D. Martin, Phys. Rev. D31, 1051 (1985). 

[12] A. G. Clark, et al, Nucl. Phys. B142, 29 (1978). 

[13] R806: C. Kourkoumelis, et al, Phys. Lett. B81, 405 (1979). 

[14] WA11: Y. Lemoigne, et al, Phys. Lett. B113, 509 (1982). 

[15] E673: S. R. Hahn, et al, Phys. Rev. D30, 671 (1984); D. A. Bauer, et al, Phys. 
Rev. Lett. 54, 753 (1985). 

[16] F. Binon, et al, Nucl. Phys. B239, 311 (1984). 

[17] E705: L. Antoniazzi, et al, Phys. Rev. Lett. 70, 383 (1993). 

16 



[18] E672: A. Zieminski, et al., Proceedings of the XXVI International Conference on 
High Energy Physics, Dallas, Texas, 1992, AIP Conference Proceedings No. 272, 
Ed. by J. R. Sanford, p. 1062. 

[19] NA3: J. Badier, et al, Z. Phys. C20, 101 (1983). 

[20] C. Biino, et al, Phys. Rev. Lett. 58, 2523 (1987). 

[21] E537: C. Akerlof, et al, Phys. Rev. D48, 5067 (1993). 

[22] B. L. Ioffe, Phys. Rev. Lett. 39, 1589 (1977). 

[23] E. L. Berger and D. Jones, Phys. Rev. D23, 1521 (1981). 

[24] E. N. Argyres and C. S. Lam, Phys. Rev. D26, 114 (1982). 

[25] E. N. Argyres and C. S. Lam, Nucl. Phys. B234, 26 (1984). 

[26] R. Baier and R. Riickl, Nucl. Phys. B201, 1 (1982). 

[27] R. Baier and R. Riickl, Nucl. Phys. B218, 289 (1983). 

[28] J. G. Korner, J. Cleymans, M. Kuroda, and G. J. Gounaris, Nucl. Phys. B204, 
6 (1982). 

[29] E. L. Berger and S. J. Brodsky, Phys. Rev. Lett. 42, 940 (1979). 

[30] A. Brandenburg, S. J. Brodsky, V. V. Khoze, and D. Miiller, Phys. Rev. Lett. 
73, 939 (1994). 

[31] K. J. Eskola, P. Hoyer, M. Vanttinen, and R. Vogt, Phys. Lett. B333, 526 (1994). 
[32] J. G. Heinrich, et al, Phys. Rev. D44, 1909 (1991). 

[33] P. Hoyer, M. Vanttinen, and U. Sukhatme, Phys. Lett. B246, 217 (1990). 
[34] E772: D. M. Aide, et al, Phys. Rev. Lett. 64, 2479 (1990). 
[35] E537: S. Katsanevas, et al, Phys. Rev. Lett. 60, 2121 (1988). 



17 



[36] E772: D. M. Aide, et al, Phys. Rev. Lett. 66, 133 (1991). 

[37] R. Gastmans and T. T. Wu, Phys. Lett. B184, 257 (1987); Nucl. Phys. B291, 
731 (1987); The Ubiquitous Photon: Helicity Method for QED and QCD, Claren- 
don Press, Oxford, 1990. 

[38] J. F. Owens, Phys. Rev. D30, 943 (1984). 

[39] J. F. Owens, Phys. Lett. B266, 126 (1991). 

[40] K. Hagiwara, A. D. Martin, and A. W. Peacock, Z. Phys. C33, 135 (1986) 

[41] J. Badier, et al, Z. Phys. C18, 281 (1983). 

[42] J. S. Conway, et al, Phys. Rev. D39, 92 (1989). 

[43] M. A. Doncheski and R .W. Robinett, Pennsylvania State University preprint 
PSU/TH/140 (January 1994). 

[44] N. Isgur and M. B. Wise, Phys. Lett. B232, 113 (1989); Phys. Lett. B237, 527 
(1990). 

[45] W. M. Tanenbaum, et al, Phys. Rev. D17, 1731 (1978); M. Oreglia, et al., Phys. 
Rev. D25, 2259 (1982). 

[46] B. Guberina, et al, Nucl. Phys. B174, 317 (1980). 

[47] P. Chiappetta and M. Le Bellac, Z. Phys. C32, 521 (1986). 

[48] M. B. Einhorn and S. D. Ellis, Phys. Rev. D12, 2007 (1975). 

[49] E. Braaten, M. A. Doncheski, S. Fleming, and M. L. Mangano, Fermilab-PUB- 
94-135-T. 

[50] D. P. Roy and K. Sridhar, CERN-TH.7329/94 (June 1994). 
[51] A. R. Clark, et al, Phys. Rev. Lett. 45, 2092 (1980) 
[52] J. J. Aubert, et al, Nucl. Phys. B213, 1 (1983) 

18 



[53] M. D. Sokoloff, et al, Phys. Rev. Lett. 57, 3003 (1986). 
[54] M. G. Ryskin, Z. Phys. C57, 89 (1993). 



FIGURE CAPTION 



Figure 1. Leading- twist predictions of the parameter A in the decay angular 
distribution of J ftp's produced in pion-nucleon collisions at = 300 GeV, plotted 
as a function of xp. (a) The three solid curves show the decay distributions of J/^'s 
produced via radiative decays of the \2 and xi states and "directly" in gluon fusion. 
The dashed curves show the effect of smearing the transverse momentum distribution 
of the beam parton by a Gaussian function exp [— (A;_i_/500 MeV) 2 ]. (b) The combined 
decay distribution of all J/V's, including contributions from xi,2 decays and direct 
production, is shown here. The lower curve shows the effect of adjusting the relative 
normalization of the different contributions to their measured values (see Table 2) by 
appropriate i^-factors. The dashed curve shows the effect of transverse momentum 
smearing and i^-factors adjustments . The data is from the Chicago-Iowa-Princeton 



(252-GeV vrW collisions, Ref. full circles) and E537 (125-GeV vrW collisions, 

Ref. JO]]; open circles) experiments. 
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